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We present measurements on parametrically driven surface waves (Faraday waves) per- 
formed in the vicinity of a bi-critical point in parameter space, where modes with har- 
monic and subharmonic time dependence interact. The primary patterns are squares in 
the subharmonic and hexagons in the harmonic regime. If the primary instability is har- 
monic we observe a hysteretic secondary transition from hexagons to squares without 
a perceptible variation of the fundamental wavelength. The transition is understood in 
terms of a set of coupled Landau equations and related to other canonical examples of 
phase transitions in nonlinear dissipative systems. Moreover, the subharmonic-harmonic 
mode competition gives rise to a variety of new superlattice states. These structures are 
interpreted as mediator modes involved in the transition between patterns of fourfold 
and sixfold rotational symmetry. 



1. Introduction 

The Faraday experiment has nowadays become a model system for pattern formation 
in hydrodynamic systems (Faraday (1831)). For a review see Miles (1990)| and Miillcr 



Friedrich fc Papathanassiou (1998)| . Standing waves are generated on the liquid air inter- 
face in response to a time periodic gravity modulation. Under typical laboratory condi- 
tions and assuming that the excitation acceleration is sinusoidal with g(t) = gp + a sin Qt 



these surface waves oscillate with twice the period of the external drive (Benjamin & 
Ursell (1954)| ). This is a consequence of the parametric drive mechanism and denoted 
here as the subharmonic Faraday resonance. Surface waves synchronous (harmonic) with 
the drive can be generated, too. They have been observed first by adding a second 
frequency component to the excitation signal ( Edwards fc Fauvc (1994)| ). Later on, fol- 



lowing a suggestion of Kumar (1996)|, harmonic Faraday waves have also been excited 
with the usual single frequency drive ( [Mullcr ct al. (1997) ). This, however, requires rather 



extreme (parameter) conditions, namely thin fluid layers in combination with drive fre- 
quencies lower than some threshold fb- Increasing / = fi/(27r) beyond /& lets the Faraday 
waves resonate with their usual subharmonic time dependence (Wagner, Mullcr & Knorr 
(2000)). For operating frequencies / close to the bi-critical value ft, the harmonic and 
subharmonic modes compete. Owing to the dispersion of surface waves different fre- 
quencies imply different wavelengths. As a consequence non-linear pattern formation is 
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Pattern formation in the Faraday instability 
Figure 1. Sketch of the experimental set-up. See text for further explanation. 



affected in a significant manner: Subharmonic modes (/ < /&) form square patterns, har- 
monic modes (/ > /&) hexagons. The occurrence of a primary hexagonal surface tiling 
of the subharmonic regime is generically due to a three wave interaction. At elevated 
drive amplitude we observe a transition towards a square pattern. This is similar to the 
canonical hexagon-line transition in Rayleigh-Benard convection, which can be observed 



if " non-Boussinesq" effects become significant (Walden & Ahlers (1981) and Ciliberto 



Pampaloni & Percz-Garzia (1988)). A transition from hexagons to squares has been found 
only recently in the Bcnard Marangoni instability ( Nitschkc fe Thess (1995)] , [Bcstchorn 
(1996)| and|Eckcrt, Bcstchorn fc Thess (1998)|). 



The measurements presented here give a comprehensive account of our investigations 
on Faraday pattern selection in the vicinity of the bi-critical point. Thereby the interac- 
tion between harmonic and subharmonic modes of different wavelengths gives rise to new 
resonant phenomena: superlattices with either fourfold or sixfold rotational invariance. 
Though superlattices are very common in solid state and surface physics, they have been 
found on macroscopic scales only recently (Pampaloni et al. (1997), [Kudrolli, Pier fe 



Gollub (1998)|, [Arbell fc Fincberg (1998)| , |Wagner, Miiller fc Knorr (1999) 
Miiller fc Knorr (2000)|) 



Wagner, 



Within a cascade of secondary phase transitions superlattices are found to mediate 
between the two incompatible symmetry classes, of squares and hexagons. For instance a 
primary subharmonic pattern with quadratic surface tiling experiences a crossover to a 



hexagonal superlattice via two quadratic superlattices (Wagner, Miiller & Knorr (2000)) 
with a prominent displacive character in one or two lateral directions. After passing 
a phase with a hexagonal superlattice the transition process reaches a pure hexagonal 
pattern characterized by a single wavelength and oscillating in synchronous response to 
the external drive. 

For several of the observed transitions we are able to provide explanations in terms 
of resonant amplitude equations for the governing spatial modes. The structure of these 
equations is simply based on symmetry and resonance arguments. In spite of their sim- 
plicity these equations provide an understanding of many remarkable features of the 
superlattices, in particular their displacive character. This phenomcnological approach 
is certainly facilitated by the small number of experimental control parameters. This is 



un like earlier experiments (Kudrolli, Pier fc Gollub (1998) and Arbell fc Fincberg (1998) 



or Wagner, Miiller fc Knorr (1999) ), which use a more complicated multiple frequency 
drive or a viscoelastic fluid to drive the system into the bi-critical situation. Thereby dif- 
ferent kinds of superlattices have been reported as well. But clearly, a larger number of 
control parameters renders a theoretical understandin g more unwieldy and les s int uitive. 
For the theoretical approach to superlattices see e.g. Dionnc & Silbcr (1997) and Silbci 
fc Skeldon (1999)| 



2. Experimental set-up 

2.1. Vibration system and sample fluid 

Figure [0 shows a schematic diagram of the experimental setup. Its heart is a large dis- 
placement shaker unit (V617 Gearing & Watson Electronics Ltd.) connected to a 4kW 
power amplifier. The shaker supplies a maximum force of 4670 N and a peak-to-peak 
elevation of s m ax = 5Amm. Such a large displacement is necessary to obtain a sufficient 
acceleration a at lower driving frequencies. The drive signal for the power amplifier is 
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synthesized by means of a digital-analog card installed in a Pentium PC. The actual ac- 
celeration of the container is measured with a piezo-electric device, the amplified signal 
of which is routed to the PC for data acquisition. Since the characteristics of the shaker 
turned out to be rather non-linear at operation frequencies below / = fi/(27r) < 10Hz 
a continuous control of the excitation signal was necessary. To guarantee a sinusoidal 
container acceleration a sin Qt the recorded accelerometer signal was decomposed into 
Fourier components. The parasitic higher harmonics of O were eliminated by admixing 
Fourier contributions with appropriate inverse phases to the excitation signal. Their am- 
plitudes were determined by a proportional control loop. That way the power spectrum 
of the accelerometer signal is made monochromatic with a purity of 99%. 

The cylindrical container for the sample liquid was machined out of aluminum and was 
anodized black. To avoid pollution and temperature drifts within the fluid, the container 
was sealed with a glass plate. The inner container diameter was d = 290mm, the depth 
50mm. Over a distance of 12mm from the edges of the container the depth continuously 
increases from zero to the bottom. This " soft boundary condition" with an average angle 
of 30° helps to minimize the generation of parasitic meniscus waves. A meniscus under 
vertical vibration always emits waves with the frequency / of the external drive. Since 
these waves have non vanishing amplitudes even at subcritical drive amplitudes a < a c 
they blur the onset detection. The beach like boundary fullfi.il their purpose well, at least 
above 10Hz. 

The probe fluid was a low viscosity Silicon oil (Dow Corning 200) with the manufacturer 
specifications of kinematic viscosity v = 5x 10~ 6 m 2 /s, surface tension a — 0. 01947V/ m 
and density p — 920kg /m 3 at our working temperature T = 25° C. A heating foil was 
mounted on the outside of the container. By means of a temperature controller the tem- 
perature measured by a PT-100 resistor (embedded in the container body) was regulated 
by ±0.1°C. 

2.2. Visualization technique 

To visualize the surface profile we used a full frame CCD camera (Hitachi KPF-1) situated 
above the fluid surface in the center of a ring consisting of 120 LEDs. The ring had a radius 
of R = 0.3m and its distance from the fluid surface was L = 1.50m. The camera was 
synchronized to the excitation signal with an exposure time of 1/256 of the drive period. 
It follows from geometrical optics that only surface elements with a certain steepness 
reflect light into the camera. 

For an evaluation of the spatial symmetry of the surface deformation £(x, y) we relied 
on a Fourier technique. To that end the recorded light intensity I(x, y) of a video image 
was convoluted with a Gaussian window function and processed by a FFT algorithm. This 
yields the two dimensional spatial power spectrum P(k). To determine the wavelength 
of the pattern P(k) is azimuthally averaged by integrating over circles with constant 
radius |k| = k. The primary peak in the resulting one-dimensional spectrum is fitted by 
a Gauss function the center of which determines the fundamental wave number. Clearly, 
the resolution of this procedure is limited by the number of wavelengths in the container. 
This is especially the case for subharmonic Faraday waves where the uncertainty of Ak/k 
is about 10%. 

Due to the nonlinear relationship between the surface elevation £(x, y) and the recorded 
light intensity I(x,y), the power spectrum entails higher harmonics of the fundamental 
wave number, even if the surface profile £(x, y) does not. Thus the relation between I(x, y) 
and 77(21, y) is generally too complicated to allow a reconstruction of the surface profile. 
Nevertheless for simple surface patterns (such as squares) we have solved this "inverse 
problem" by the following method: Starting from an estimated surface profile composed 
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Figure 2. Neutral stability curves a(k) computed for the parameters of the sample fluid at a 
drive frequency a)/ = 6.25Hz < fo and b)/ = 7.25Hz > /;,. In a) the primary resonance is 
harmonic, in b) it is subharmonic. Regions where the flat surface state is unstable are shaded. 
Horizontal lines denote the thresholds for secondary and higher order transitions (for details see 
text). 



of a small number of spatial Fourier modes, the light distribution of the expected video 
image was computed by means of a ray tracing algorithm. Then we adapted the mode 
amplitudes and their relative phases such as to optimize the agreement between the 
calculated and recorded video picture. 

A reconstruction of the full time dependence of an oscillating surface wave pattern 
was not possible with our equipment. Nevertheless, the electronic shutter of the camera 
provides an easy and very sensitive technique to discriminate subharmonic frequency 
components in an otherwise harmonic time signal. This is because a harmonic time 
dependence £ (£) is invariant under the symmetry operation t — ► t + ^ implying a 
frequency spectrum of integer multiples of O thus ( h (t) — Yl n Cn& mQt ■ In contrast the 
subharmonic time signal transforms after one drive period as C s ~~ * — C s enforcing a 
Fourier representation in the form of £ s (£) = X) n Cn e ■ Thus by triggering 

the camera shutter with the drive frequency Q, video images with a harmonic time 
dependence appear stationary, while those with subharmonic frequency contributions 
flicker due to a slight optical asymmetry between heaps and hollow of the deformed 
surface. Note however, that this trigger technique does not allow to identify harmonic 
frequency components in an otherwise subharmonic spectrum. 



3. The onset of the Faraday instability 

It is well known that the stability problem of a free liquid surface under gravity mod- 
ulation (Faraday i nstability) can b e approximately mapped to that of a parametrically 



driven pendulum ( |Raylcigh (1883)| , |Bcnjamin fc Urscll (1954)| ) . The primary resonance 



of which occurs at twice the per iod of the drive (subharmonic response). However, as first 



pointed out by Kumar (1996) the Faraday instability may also appear in synchronous 
resonance with the external drive, usually denoted as the harmonic response. The con- 
ditions under which the harmonic resonance preempts the subharmonic one have been 



worked out in detail by Cerda fc Tirapegui (1997) and Muller et al. (1997) revealing that 



low filling levels in combination with small drive frequencies are necessary. In the present 
experiment we choose a fill height of ft. = 0.7mm, which is - at the operation frequencies 
of 6 < / < 8Hz - comparable to the viscous penetration depth £ = \f2vjQ, m 0.5mm. 
For the fluid parameters at hand a linear stability analysis of the flat surface state (ac- 



cording to the method of Kumar & Tuckcrman (1994), which assumes a laterally infinite 



system) reveals the location of the bi-critical point at a drive frequency of /;, = 6.3Hz. 
Figure @ shows neutral stability diagrams (drive amplitude a vs. wave number k) for 
both situations f < ft and f > ft- Experimentally the critical acceleration a c (absolute 
minimum of the neutral stability diagram) has been determined by setting up the system 
at a constant frequency and ramping a quasi-statically in steps of 0.2% suspended by 
intervals of 240s. The onset amplitude a c was defined when the camera detected the first 
light reflex (figure Cm). To enhance the detection sensitivity the surface was illuminated 
by a diffusive light source from the side rather than using the dark-field technique de- 
scribed above. This is because the latter method requires a minimum surface gradient 
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Figure 3. (a) Critical amplitude a c and (b) critical wave number k c for the onset of the Faraday 
instability drawn as a function of the drive frequency /. The bi-critical point / = ft is located 
where the harmonic (/ < ft) and subharmonic (/ > /&) thresholds intersect. Symbols mark 
experimental data points, lines the theoretical results for a laterally infinite fluid layer. Circles 
and dotted lines refer to the harmonic response, squares and solid lines to the subharmonic one. 



Figure 4. Phase diagram of the observed patterns obtained by quasi-statically ramping the 
driving force. The symbols mark the observed transition points between the different patterns, 
the lines are guides for the eyes. The spatial ordering of the patterns is indicated by Roman 
numerals, the arabic letters s and h denote the character of the time dependence being either 
purely subharmonic or harmonic. With s+h we indicate patterns formed by a interaction of 
subharmonic and harmonic modes. The thick line separates harmonic from subharmonic and 
s+h regions. I: flat surface; II: subharmonically oscillating squares (figure nil a,b) ; Ilia: V2x y/2 
superlattice p2mg (figure hlc,d); Illb: y/2 x \/2 superlattice c2mm (figure hlkf); within the 
subregion above the dotted ine the pattern is time dependent and disordered (see figure hj) ; 
IV: \/3 x \/3 superlattice, (only for / < 6.9Hz stationary, figure [la): V: harmonically oscillating 
hexagons (figure Ha); VI: harmonically oscillating squares (figure Bid); VII: 2x2 superlattice 
(figure hq); VIII: local instability and droplet ejection. 



of |V£(ir,y)| = tana w 0.1 for the onset detection. We estimate the accuracy of our 
threshold determination by 0.5%. 

Once a standing wave pattern had covered the whole surface the fundamental frequency 
of the surface oscillation was determined with the help of the electronic shutter of the 
video camera. That way we located the transition point at a bi-critical frequency fb = 
6.5 ± 0.1Hz. After these preliminary measurements we switched back to the dark-field 
illumination to proceed with the spatial pattern analysis. The critical wave numbers 
fc^ and respectively k s c were determined by Fourier transforming a surface image taken 
at a driving strength of e = (a — a c /a c ) « 3% (figure |3Jb). The operating prescription 
k c = k(e ~ 3%) for the determination of the critical wave number is motivated by the 
fact that we were unable to detect any change of the wave number k by varying e (sec 



also Wernet et al.(2001) ). The experimental results for a c and k c as well the bi-critical 
frequency fb are (found to be) in good agreement with the theoretical predictions. For 
the critical acceleration the discrepancy is less than 2%. Here the uncertainty is mainly 
due to errors in the determination of the small fill height h. For larger values of h the 
agreement improves up to 1%. For the critical wave number the discrepancy between 
theory and experiment is better than 4% in the harmonic case but it increases up to 10% 
on the subharmonic regime. This is due to the spatial resolution, which becomes worse 
at larger wavelength. Owing to the abrupt change of the response frequency at / = fb 
the wave number shows a discontinuous jump (see figure pp). The empiric ratio of the 
wave numbers at / = fb is found to be 

^-U P = 1.58 ±0.15, (3.1) 

k h 

in agreement with the prediction of the linear stability theory p-|th eo = 1.59. 



4. Overview over the phase diagram 

The phase diagram shown in figure has been obtained at various constant driving 
frequencies / = f2/(27r) while ramping the driving amplitude from e = —1% up to 25% 
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in steps of 0.2%. After each increment the ramp was suspended for 240 seconds to give 
the system time to relax. Then a photo or - in the case of time dependent patterns a 
video film - of the surface state was taken. At the point where a new spatial or temporal 
mode appeared or an existing one died out, the actual acceleration was defined as the 
transition boundary to a new " phase" . At the maximum acceleration amplitude the ramp 
was reversed to check for an eventual hysteresis. 

In the following chapter 5 we describe in detail the type of primary patterns, which 
appear near onset of the Faraday instability. Secondary and higher order transitions 
towards more complicated structures are dealt with in sections 6 and 7. Thereby two 
representative experimental runs will be described in detail, the first was taken at / = 
6.25 Hz< fb and a second at / = 7.25 Hz> /&. 

In the former case the primary pattern exhibits a harmonic time dependence, which 
turns out to be quite robust as it persists over the whole investigated e-ramp. The primary 
spatial surface wave structure starts with an ideal hexagonal symmetry (region V) , which 
then transforms into a pattern of squares (region VI) as e is raised. This transition is 
hysteretic, its global aspects can be understood in terms of a simple model of six coupled 
amplitude equations. 

In the second run at f > fb the primary surface pattern consists of subharmonically 
oscillating squares (region II). On increasing the drive strength e the interaction with the 
neighboring harmonic Faraday instability leads to the appearance of a quadratic v2 x V2 
superlattice with a displacive character in one and/or two lateral directions (regions Ilia 
and Illb, respectively). After crossing a phase region of non-stationary patterns with a 
slow time dependence the system enters a hexagonal V3 x -\/3 superlattice (region IV). 
Mediated by a second local reconstruction process the final stationary surface pattern is 
a quadratic 2x2 superlattice (region VII) . Regardless of whether the response is s or h, 
the surface finally breaks up and droplets are ejected (region VIII, e ~ 15 — 25%). 



5. Pattern formation close above onset of the Faraday instability 

Close to the harmonic onset of the instability (/ < fb) hexagons are the preferred 
primary surface pattern (p|a, regio n V) for f > fb, however, squares are stable (figure |ll|a). 



Wagner, Miiller fe Knorr (2000) (see figure 3 in it) has shown that even for small e the 



wave profile is rather anharmonic. 

5.1. Theoretical model for the primary hexagonal pattern at f < fb 

The appearance of hexagons at the " harmonic side" of the bi-critical point follows from 
a triple wave vector resonance: Since the spectrum of the harmonic Faraday mode C h (i) 
consists of integral multiples of the drive frequency it allows for two critical standing wave 
modes ki and k.2 to resonante with a third one. The requirement |ki| = Jlcgj = |lt3 1 = k'* 
along with the resonance condition ki + k2 + k3 = enforce a mutual angle of 120° 
between the wave vectors implying the hexagonal symmetry. The evolution equations for 
the respective mode amplitudes H± , H2 , H3 are of the following structure 

8 t H 1 =eH 1 +0H$H$- [{H^ 2 +T(120°)(\H 2 \ 2 + \H 3 \ 2 )] H,. (5.1) 

Thereby j3 is a second order coupling coefficient and Y{&) is the cubic cross coupling 
coefficient, which depends on the angle between the interacting modes. Moreover the 
stars denote complex conjugation. The corresponding equations for H\ and H2 follow by 
permutation of the indices. The term of cubic order is crucial for saturation. A linear 
stability analysis of the finite amplitude stationary solution \H±\ = \H.%\ = \H 3 \ given 
by |Cifibcrto, Pampaloni fc Pcrez-Garzia (1988) yields a backwards bifurcation out of 
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Figure 5. Photographs of the fluid surface as obtained by ramping the drive amplitude e at 
/ = 6.25Hz. In the corners of the pictures the container boundary is visible. In (a) and (d) 
circles mark the nuclei at which the transition process is initiated, (e) is an enlarged sector of 
(b) with marked penta-lines and hepta-defects. 

the trivial solution Hi = 0. This reflects a hysteretic transition from the undisturbed 
flat surface to a pattern of hexagons (compare figure M). However, we were unable to 
resolve any hysteresis because of small amplitude (harmonically oscillating) meniscus 
waves emitted from the rim of the container. 

5.2. Theoretical model for the primary square pattern at f > ft, 
Understanding the pattern selection process at the subharmonic side of the bi-critical 
point is rather more complicated. Since the frequency spectrum of the subharmonic 
Faraday response consist of half integer multiples of /, any triple of linear unstable modes 
is prevented from resonating. Thus nonlinear pattern selection is dominated by a four- 
mode interaction. Unlike triad resonances, which operate exclusively at an interaction 
angle of 9 = 120°, four- wave resonances are less selective as they work at arbitrary angles 
8. This fact is also reflected by the corresponding system of amplitude equations. Taking 
a set of N standing waves with wave numbers k^ at length |kj| = kl but arbitrary relative 
orientation, the respective mode amplitudes S% are governed by the following evolution 
equations 

N 

d t S i =eSi-J2 r (0*3)\ s i\ 2 Si, ( 5 - 2 ) 

.7=1 

with 9ij being the angle between ki, k,. Usually the participating modes are taken to 
be equidistant on the circle |kj| = fc*, thus 0j,i+i = 2ir/N. In this case iV indicates the 
type of symmetry of the pattern, namel y N = 1 lines, N = 2 squares , N — 3 triangles 
or hex agon s, ... . As outlined in Refs. ( Milner (1991)| , Muller (1994) . Zhang fe Vinals 



(1997) and Chen & Vinals (1997)) the question of what is the most preferred symmetry 
is reduced to minimizing the "free energy" 

N N 

F = -e^Sil + 2 E r Vii) l^l 2 l^l 2 ( 5 - 3 ) 

with respect to N at given T(8). For low viscosity fluid layers of infinite depth the coupling 



function T(9) was first evaluated by Zhang fc Vinals (1997) , who found that a pattern 



of square symmetry is the most preferred one at drive frequencies beyond / ~ 50Hz. 
At lower frequencies patterns with a degree of rotational symmetry up N = 7 (quasi- 
periodic) are likely to occur These predictions were found to be in qualitative agreement 



with experiments (Kudrolli fc Gollub (1996), Binks fc van de Water (1997)). The latter 



authors also extended the above considerations to the case of finite fill heights and found 
square patterns to dominate also at lower drive frequencies in agreement with their own 
experiments and also with ours. 



6. Secondary and higher order transitions at / < /& 

6.1. Hexagon-square transition 

In this section we investigate the crossover from the primary hexagonal structure (region 
V in figure 0) to the square pattern of region VI of the harmonic regime. Throughout 
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Figure 6. (a) Surface area covered with squares Ai related to surface area of squares and 
hexagons As as function of the drive strength e at the frequency / = 6.25 Hz. (b) The auto- 
correlation function (7(0 = 90°) of the power spectra P(k) on the circle (jfcj = k c ). Upright 
(reversed) triangles refer to the upwards (downwards) amplitude ramp. 



Figure 7. Azimuthal autocorrelation C(<f>) derived from the power spectra P(k) at the wave 
number |k| = k c as a function of the drive amplitude e for ramping up and down. 



the whole bifurcation cascade the time dependence is purely harmonic without percepti- 
ble subharmonic frequency contributions. The results described below were obtained by 
ramping the drive amplitude e = —0.1% up to e — 16.5% while keeping the frequency 
/ = 6.25-ffz fixed (see figure ||). The transition is connected with a strong hysteresis 
giving an overlap between the two ideal structures. Since the lateral aspect ratio (con- 
tainer size over wavelength) is not large the hexagonal pattern adapts to an azimuthal 
symmetry. This affects the defect dynamic of the pattern and thereby also the transition 
process. Starting the e-ramp at the hexagonal structure, a rearrangement of the pat- 
tern sets in, initiated by six nuclei of local quadratic order as indicated by the circles 
in figure Ha. The size of these patches increases with e (Fig. 0b). Conversely, starting at 
elevated e with the perfect square pattern and ramping downwards generates four nuclei 
of local hexagonal order (Fig. ||d). Along the domain boundaries between the patches 
pcnta-hepta defects occur. This kind of defect is very common in 2-D hexagonal pat- 



terns. Experimental (Ciliberto, Pampaloni & Perez-Garzia (1988) and Fam, Ohata & 
Wu (2000)) and theoretical (Prismcn & Ncpomnyashchy (1993) and Tsimring (199G)|) 



investigations reveal that a penta hepta defect can be formed by a phase defect among 
two of the participating modes. Figure |5je presents an enlarged subrange of figure |5jb 
depicting "penta lines". This is a row of penta defects (unit cells with five neighbors) 
ending in hepta defects (unit cells with seven neighbors). Strengthening the drive makes 
the domain walls invade the areas of hexagonal symmetry with new quadratic cells being 
generated along the penta lines. 

To quantify the hysteresis of the hexagonal quadratic reconstruction we applied two 
different techniques, one in real space and the other in Fourier space. In the former case 
the number of unit cells with four neighbors A± and those with six neighbors Aq was 
counted. The result of this procedure is shown in figure pk. It reveals a hysteresis loop 
extending from e = 5%± 1% up to e = 15% ± 1%. The obvious staircase behavior reflects 
the discretization of the k- values due to the finite size of the container. For different runs 
the steps do not occur at the same e-values. 

The second method for evaluating the square- vs. -hexagon surface coverage we Fourier 
transformed the full surface picture and evaluated the power spectrum P(k) = P(k, <j>) 
at |k| = k c as a function of the azimuthal angle tp. We then computed the azimuthal auto 
correlation C(</>) as follows 

cm = wv^mv ■ (6 - 1) 

That way a pattern with quadratic symmetry leads to a peak at <fi = 90° and 180°, while 
a hexagonal pattern produces maxima at <f> = 60°, 120° and 180°. Figure H illustrates 
how C((j)) develops with the drive amplitude, both at increasing and decreasing e. In 
order to compare with the cell counting method (see figure pk) the value C(</> = 90°) is 
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Figure 8. Set of wave vectors necessary to build up patterns of rolls, hexagons or squares. 



Figure 9. Sketch of the bifurcation diagram for the transition from a hexagonal to a quadratic 
pattern. Dotted lines indicate unstable, straight lines stable branches. The hysteresis observed in 
the experiment may be attributed to the bi-stable region, where squares and hexagons co-exist. 

plotted versus e in figure gb. Regarding the width of the hysteresis loop, the two methods 
agree within a few percent with each other. 

An important feature of the above phase transition is the constancy of the wave number 
k h (e) = fc'': Within the experimental resolution of Ak/k = ±1% no dependency of the 
wave number on e could be detected throughout the whole investigated drive amplitude 
range. This is of particular significance as it allows for to describe the global aspects of 
the phase transition in terms of a simple model as given in the next section. 

6.2. Comparison with theory 
It was outlined in § 5 that the hexagonal symmetry at small e is a consequence of a 



three wave resonance, reflected by the second order term in equation (5.1). However, 
upon increasing the control parameter e the term of cubic order becomes increasingly 
important thus finally enforcing the transition towards squares. A minimal model that 
allows stationary solutions in form of squares and hexagons (and also lines) along with 
a linear stability analysis of these solutions can be found in a recent publication of 
Regnier et al. (1997)| . This model is an extension of equation (pTJ) as it relies on six 
independent modes ki t i=i...6 rather than only three (see figured). Since 1^=1,2,3 enclose 
mutual angles of 120° the remaining wave vectors 1^=4 se describe a second pattern 
of hexagons, which is rotated by an angle of 30° relative to the first. Depending on 
the number of saturated mode amplitudes Hi the model allows stationary solutions in 
form of lines, squares, hexagons or even more complicated solutions without translational 
symmetries (quasi-periodic structures) . We shall disregard these complications here and 



focus on the square-hexagon competition. Regnier et al. (1997) demonstrate that the 



hexagonal state (given by finite \H±\ — \H%\ = I Hj, |) becomes unstable against shear 
distortions 6 in form of rhombuses (8H2 ^ 0^ 8H3) at a drive amplitude e > en- By 
way of contrast, for a critical value of e < es squares (given by finite |-ffi| = |H4|) become 
unstable against rhomboedric {H2, H3) disturbances. These results are summarized in 
the bifurcation diagram of figure |[ 

It is tempting to attribute the bi-stable square-hexagon region between es < e < £h 
(see figure ||) to the hysteretic region depicted in figure ||b,c). Note however, that the ob- 
served transition runs through a reconstruction via penta-hepta defects, the complicated 
space dependence of which goes beyond the scope of the present model. 

We mention that a discussion of a transition between hexagons and lines in terms 
of the above three m ode model (|5.l|) has been given earlier in the context of Raylcigh- 
Benard convection by Waldcn fc Ahlcrs (1981)| and Cilibcrto, Pampaloni & Perez-Garzia 



(1988). Thereby the phase transition results from non-Boussinesq effects induced by a 



strong applied temperature gradient. Recently, several authors report Benard-Marangoni 



experiments, which show a transition from hexagons to squares (Nitschke & Thess 



(1995),|Bestehorn (1996) and Eckert, Bestehorn & Thess (1998)). A similar transition is 



found by Abou, Wcsfreid & Roux (2000) on the surface instrability of magnetic liquids 
(Rosensweig-Instability) . However, unlike our measurements, all of the afore mentioned 
experiments with a hexagon to square transition exhibit a considerable nonlinear wave 
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Figure 10. Phase diagram of an amplitude scan at / = 7.25 Hz. Time dependent transient 
means an uncorrelated pattern (figure hj), quasistationary indicates a pattern of an almost 
perfect \/% x \/3 superlattice (compare figure h5J) with a slow defect dynamic. 



Figure 11. A 10 x 10 cm sector of the photographs of the fluid surface and corresponding Fourier 
spectra at / = 7.25 Hz. a,b): e = 3%, region II, pattern of squares; c,d): e = 6.6%, region Ilia, 
V2 x y/2 superlattice p2mg with a displacive character in ^-direction; e,f): e — 14% \[2 x v2, 
region Illb, superlattice c2mm with a dispalzive character in x- and y-direction. The dotted 
lines and arrows on the left pictures indicate the directions in which the rows and columns of 
elevation maxima are displaced. 



Figure 12. Particle model of the a) \[2 x y/2 p2mg, b)^2 x y/2 c2mm superlattice. Also 
marked are the mirror (m) and glide (g) planes. 

number variation k(e) up to 10%, which rules out a description in terms of space inde- 
pendent amplitude equations of the type given above. 



7. Secondary and higher transitions at / > /& 

In the same way as in the preceding paragraph we now turn to the bifurcation scenario 
at the opposite side of the bi-critical point at / > /&. An amplitude ramp taken at 
/ = 7.25Hz serves as a representative example. § 7.4 is an exception, where the focus is 
on the frequency regime 6.6Hz < f < 6.9Hz. 

An overview over the transition behaviour at / = 7.25Hz is given by figures [L0| and 
|ll| . Starting from subharmonic oscillating squares (region II, figure |ll|a) a transition to 
a quadratic V2 x v2 superlattice with a displacement of neighbouring elevation maxima 
in the lateral x-direction (region Ilia, figure Oc) takes place. The next pattern is again a 
V2 x v2 superlattice of the original square lattice but this time it exhibits a displacement 
in both x and y direction, (region Illb, figure Oe). After a time dependent transient (see 
figure [141 for a snapshot) this pattern transforms into a "quasi-stationary" hexagonal 
\/3 x V3 superlattice (region IV, figure 03) . With " quasi-stationary" we indicate that 
the pattern is slightly disturbed by defects, which induce a slow time dependence on the 
scale of minutes. Further raising the drive amplitude e at / = 7.25Hz makes the quadratic 



symmetry reappear in form of a 2 x 2 superlattice (region VII, figure 16). Performing 
the same amplitude ramp at a lower frequency of 6.6Hz < / < 6.9Hz the V3 x %/3 
superlattice directly reduces to its underlying pure hexagonal tiling. This last transition 
will be discussed in § 7.4. 

7.1. The quadratic \/2 x y2 superlattice 
The bifurcation sequence starts at the primary ideal pattern of subharmonically oscillat- 
ing squares (region II and figure Oa,b.) composed of the two fundamental wave vectors 
ksi and ks2- Increasing the drive strength e beyond 5% displaces every other column 
of elevation maxima in the direction ±(ksi + ks2)- This is shown by the arrows in fig- 
ure [ilia. The resulting pattern (phase region Ilia) is depicted in Oc. The displacement 
is accompanied by the simultaneous appearance of the modes km, k#2, and km in the 
power spectrum (see figure |Tl|d). Due to the approximate equality |k^i| = |k# 2 | ~ /c'' 
these modes are only slightly damped. In contrast, |krji| is significantly smaller than 
fc^and k^. Therefore the mode D\ is strongly damped. Its mission is to act as a mediator 



10 C. Wagner, H.-W. Mutter and K. Knorr 

mode, enabling the resonance between S- and iJ-modes according to the geometrical 
rules k/fj = (ks 1 + kcj and k# 2 = (kg 2 + koj). This implies that kni must be associ- 
ated with a subharmonic time dependence. Moreover, it follows from k^:/k s c ~ 1.58 that 
|k_Di| = 0.71 ~ l/v2) which leads to the commensurate relationship k^i ~ ^(kgj +kg 2 ) 
for the new fundamental wave vector. The corresponding elementary cell (see figure 13a) 
is rotated by 45° relative to the original quadratic grid, and the basic wavelength is 
enlarged by a factor of \[2. We therefore denote this pattern as a \/2 x \f2 superlattice 
(more precisely: (\/2 x \/2)R45 or \[2 x \[2 p2mg in the nomenclature of space group 
theory). It turns out that the displazive character of the superlattice is inherent in the 
phase information carried by the participating modes. Barring higher spatial harmonics 
the space dependence of the surface deformation shown in figure |ll| can be expressed as 

2 2 

n(r) = y^ Si cos (kg^r + ifej +Dicos(k Dl r + <AD 1 ) + ^ Hi cos (k Ht r + <p Hi ) . (7.1) 

i=l t=l 

By means of the ray tracing technique outlined in § 2.2 we can now simulate the video 
image associated with rj(r) and adapt it to the empiric result. Taking (f>s 1 — 4>s 2 — (by 
a proper choice of the origin) our investigation reveals that the displacement visible in 
figure Ob can only be reproduced if the spatial phases (f>r> 1 , <$>h x ? an d 4>h 2 adopt values 
close to = 7r/2. The associated surface pattern exhibits a 180° rotational symmetry. 

By increasing the drive further the y2 x V2 p2mg superlattice undergoes a transition 
which restores the fourfold symmetry. Similar to the above described shift of the columns 
of elevation maxima, it is now additionally the rows, which experience a displacement 
in the direction ±(ks 1 — ks 2 ) (indicated by the arrows in figure |ll|c). In figure O we 
measured the amount of symmetry restoration by comparing the spectral power associ- 
ated with D\ and D 2 - Beyond e ~ 9% the transition is complete. The resulting so called 
y/2 x \/2 c2mm superlattice (see figure |l2|b) is depicted in figure |ll| e and associated with 
the phase space region Illb. In Fourier space this transition is carried by the additional 



modes k# 3 , k# 4 and k£> 2 as shown in figure |ll| f. Extending the surface representation 7T 
by these additional components and using it to re-perform the ray tracing image analysis 
yields the phase information <j)£ )i i=1 2 = (f>H t i=1 4 = it/2. 

We mention that the y/2 x \J~2~ c2mm superlattice state was not observed in an earlier 



measurement on a liquid of higher viscosity (Wagner, Miiller & Knorr (2000)). Thereby 



this structure is preempted by a transition to a hexagonal symmetry. 

7.2. Theoretical model for the displacive phase 

The experimental investigations outlined in the previous section reveal that the promi- 
nent displacive character of the y/2 x y2 superlattices is associated with the phase infor- 
mation carried by the spatial Fourier modes. In what follows, a minimal model is con- 
structed, which is able to explain the experimentally measured phases. It is important 
to point out that the structure of these equations just relies on symmetry and resonance 
arguments (triad wave vector resonances). The numerical values of the appearing co- 
efficients are not known; their evaluation would require a rather complicated nonlinear 
analysis. Furthermore, for the sake of simplicity, we limit our discussion to the \/2 x y2 
p2mg superlattice. The generalization to the more symmetric y2 x y2 c2mm pattern is 
straightforward . 

Assuming that the amplitudes of the primary square pattern have settled at some finite 
value 5i = S*2 7^ the leading order behavior of the remaining modes is governed by the 
following set of equations 

d t D 1 = e D D 1 + li d {Stm + S* 2 H 2 ) + xsSiS 2 Dl 
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Figure 13. Integrated intensity of the Di-peak (circles) and the D2-peak (squares) in the power 
spectrum as a function of the drive amplitude e. The error bars mark the standard derivation 
resulting from 5 succeeding runs. 



Figure 14. Snapshot of the fluid surface as obtained at / = 7.25 Hz and e = 17%. The 
pattern is strongly time dependent and spatially uncorrelated. 



Figure 15. A 8 x 8 cm sector of the photographs of the fluid surface and corresponding Fourier 
spectra at / = 6.65 Hz and e = 5%. This photograph has been obtained by using a diffusive light 
source mounted off axis instead of the standard visualization technique. That way the character 
of the hexagonal Vo x \/i superlattice (region IV) appears to be more pronounced (compare 
figure 4e in (Wagner (2000))). 

d t H x = e H H x + fi H S 1 D 1 (7.2) 

d t H 2 = e H H 2 + n H S 2 Di. 

Thereby Eh < and £d<0 are the coefficients of linear damping of respectively the H 
and D modes, while [1d,h and xs are nonlinear coupling coefficients associated with the 
triad wave vector resonance. 

Although the model equations ( |7.3| ) are linear in D and H and thus saturation is 
not implied, the appearing nonlinearities are phase selective. By writing the complex 
amplitudes in the form Aj = |Aj| exp (i^,^) with again <pgi — 4>S2 taken to be zero 
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(7.4) 



with m being an integer. By inspection one finds that the solution with even m leads 
to an square pattern with amplitude modulation while odd m gives rise to the observed 
displacement (phase modulation). Which one is stable, depends on the numeric values of 
the coefficients. From the experimental data we conclude that <^h\ = (f>£) i = ir/2 is the 
solution applicable to the present experiment. 

7.3. The hexagonal \/3 x \/3 superlattice 

Upon further increase of the driving force at / = 7.25Hz the quadratic v2 x \[2 c2mm 
superlattice transforms into a "quasi-stationary" hexagonal superlattice (region IV, com- 
pare figure jl5|) after passing a region of transient time dependent (figure [m]) with squares, 
hexagons and disordered patterns appearing at the same value of e. The term "quasi- 
stationary" is to indicate that the pattern is affected by defects on a slow time scale of 
minutes. The power spectrum shown in Figure [lqb reveals that the structure is composed 
of a set of three harmonic modes Hj and three subharmonic modes Sj . The angles be- 
tween the wave vectors of the three H modes is 120° and so is the angle of the S-modes. 
The H-pattern is rotated by 30° and locked in phase with respect to the S-pattern, the 
ratio kh/k s is 1.73 ~ V3. Therefore the pattern is denoted as a V3 x V3 superlattice. We 
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mention an earlier investigation on a more viscous fluid of v = lOcS ([Wagner, Miiller fc 



Knorr (2000)|): Thereby the transition from the quadratic \/2 x \/2 p2mg to the hexag- 
onal \/3 x \/3 superlattice took place in a more correlated manner via stacking faults 
while the y/2 x \/2 c2mm state could not be observed. 

At higher e-values the further development of the ^/3 x v3 superlattice depends on the 
drive frequency: At 6.6Hz< / < 6.9Hz, where the superlattice is perfectly stationary and 
almost free of defects, the amplitude of the subharmonic components of the structure con- 
tinuously become smaller until the hexagonal base pattern with a pure synchronous time 
dependence remains (region V, compare figure J5|a). This trannsition is slightly hysteretic 
(Ae ~ 1 — 2%) and can be localized accurately by means of the triggering technique 
described in § 2. We come back to this crossover in § 7.4, in order to present a theoretical 
model. 

At drive frequencies larger than ss 6.9Hz (for concreteness let us return to our run at 
/ = 7.25Hz) the hexagonal y/3 x -\/3 superstructure transforms into a 2 x 2 superlattice 
(region IV — v VII) thereby restoring the quadratic symmetry. The crossover takes place 
via a spatially weakly correlated transient, being subjected to a strong temporal dynamic. 
The corresponding pattern looks similar to the surface state shown in figure |IJ. In its 
final state the subharmonic and the harmonic wave vectors \ts,H are aligned with each 
other (figure |bqb). With k h rs k^ but k s rs 0.8fc^ the resulting length ratio kn/ks is 
about 2. 

7.4. Theoretical model for the v3 x v3 superlattice-to-hexagon transition 

The following section provides a theoretical model for the transition from the V3 x v3 
superlattice (region IV) to the pure hexagon state in region V, taking place at 6.6Hz 
< / < 6.9Hz. Following the lines of § 7.2 spatial and temporal resonance arguments are 
used to construct the set of amplitude equations. We consider the underlying hexagonal 



base pattern, carried by the modes Hi (c.f. equation (5.1)) as being saturated at some 
finite amplitude. Then the bifurcation of the subharmonic modes Si, S 2 , S3 is governed 
by the following system of equations 

3 
dtSi = eSi + [3 HS (HiS 3 + H*S 2 ) - /^ |#/ 5 i + h -°- t - ( 7 - 5 ) 

i=i 
where the equations for S 2 and S3 follow by cyclic permutation, and (3hs an d /1 are 
nonlinear coupling coefficients. In this model the resonance condition between the S and 
the H modes enforces a wave number ratio kn/ksltheo — V3 « 1.73. Our measurements 
reveal \k h \ = kf (e = 4%) and |k s | = 0.91fcf giving k H /k s \ex P ~ 1-73. This justifies 
our assumption that harmonic modes dominate the pattern forming process, while the 
(slaved) subharmonic ones have to adapt their wavelength. 

The occurrence of the V3 x v^ superlattice is caught by a linear stability analysis of 
equation ( frT5| ) . To that end it is convenient to take the amplitude H = \Hi\ = \H 2 \ = | H3 1 
as the control parameter. Then one obtains a positive growth rate (instability) for the 
S t if 

e + 2(3 H s\H\-^\H\ 2 >0. (7.6) 



This condition determines the transition to the superlattice. Recall that (7Jj) does not 
entail terms which are nonlinear in Si thus the model neither explains saturation of the 
Si nor the experimentally observed hysteresis in the transition. 

Let us now focus on the spatial phase dynamic of the hexagonal superlattice. For 
the underlying hexagonal structure with amplitudes Hi = \Hi\ exp («$#;) we have ei- 
ther §h = >~2j 4*Hj = or 7T, corresponding respectively to "up" or "down" hexagons. 
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Figure 16. A 8 x 8 cm sector of the quadratic 2x2 superlattice and corresponding Fourier 
spectra at / = 7.25Hz and e = 24% (region VII). 

For Faraday waves this distinction is not significant because the surface oscillation pe- 
riodically switches between these two possible states. Splitting the subharmonic mode 
amplitudes also into modulus and phase and assuming without loss of generality that 



4>Hi = 4>n-i = 0, the imaginary part of equation (7.5) yields the identity (j>s t — 4>s 2 = 4>S 3 - 
In order to determine that value it is necessary to proceed with the amplitude expansion 
up to the quintic order. Barring all terms which are not phase selective the extension of 
equation ( |7.5| ) reads 

dtSi = sSi + (3ns (H X S 3 + H*S 2 ) + .. + xSi (Stf (S* 3 f + ..., (7.7) 

with x being a no nlinear coefficient. Introducing the parametrization $5 = J^ ■ <pg i = 



2>4>Si equation (7/7) yields the stationary solution 

$ s = rm/2. (7.8) 

where n is an an integer. Following the usual nomenclature the pattern with odd n is a 
triangular superlattice, as its rotational symmetry is threefold. For even n the superlattice 



exhibits a sixfold symmetry (c.f. Muller (1993)). Since a sixfold symmetry center can be 



easily identified in figure p_5L we conclude that n is even in our experiment. We point 
out that an example of a \/3 x v3 with third order rotational symmetry (n is odd) has 



recently been observed by [Pi et al. (2000)| (see also |Wagner fc Muller (2001)| ) 



8. Conclusions 

We have presented a comprehensive investigation on Faraday wave pattern formation in 
the vicinity of a bi-critical situation. Thereby modes with harmonic and subharmonic time 
dependence interact and lead to a variety of superlattice states. This kind of structure 
acts as a mediator during the transition process between two incompatible space groups 
(squares and hexagons). We developed several sets of amplitude equations to model the 
observed phase transitions. Special attention is devoted to the phase information carried 
by the participating modes, which is responsible for the remarkable displacive feature of 
one of the superlattices. 

Superlattices are rather common in 2-D solid state physics, and a comparison is there- 
fore instructive. The transition from a simple hexagonal lattice to a -\/3 x \/3 super- 
structure has been observed for instance in monolayers of C2CIF5 adsorbed on graphite 
( Fassbender et al. (1995)| ). The transition from the subharmonic quadratic base pattern 



to a v2x v2 superlattices with a displacive character in one direction (p2mg, region Ilia) 
i s analogous to the rec onstruction of the clean (100) surface of W (=tungsten) crystals 
( |5chmidt et al. (1992) ). Here the surface atoms are displaced in exactly the same way 



as the elevation maxima of the surface profile in the present study. Most interestingly, 
the surface of W crystals obeys a transition to a \[2 x \[2 superlattices with a displacive 
character in two directions (c2mm, identical to our pattern in region 111b) in the presence 
of Hydrogen atoms. 

We were also able to provide for the first time a pattern forming system, which under- 
goes a hexagon-to-square transition without a simultaneous change of the fundamental 
wavelength. This is in contrast to earlier observations on the Benard Marangoni system, 
which showed considerable wave number variations during the transition process. The 
accompanying description in terms of well defined spatial Fourier modes, makes the phase 
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transition very fundamental and identifies the underlying pattern selection mechanism 
to switch from a triad to a four-wave vector resonance. 

We thank J. Albers for his support. This work is supported by the Deutsche Forschungs- 
gemeinschaft. 
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